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ABSTRACT 

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral 
energy distributions of 63 bright (L(Lya) > 10 43 erg s l ) Lya emitting galaxies (LAEs) in the redshift range 
1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders 
of magnitude, from 7.5 < lo gM/M Q < 10.5. Moreover, although most LAEs have small amounts of extinction, 
some high-mass objects have stellar reddenings as large as E(B — V) ~ 0.4. Interestingly, in dusty objects the 
optical depths for Lya and the UV continuum are always similar, indicating that Lya photons are not undergoing 
many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can 
vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star-formation-rate-log-mass 
diagram, our LAEs fall above the “main- sequence” defined by z ~ 3 continuum selected star-forming galaxies. In 
this respect, they are similar to submillimeter-selected galaxies, although most LAEs have much lower mass. 

Key words : cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: starburst 
Online-only material: color figures 


1. INTRODUCTION 

Patridge & Peebles (1967) originally predicted that the Lya 
emission line could be a very useful probe of the high-redshift 
universe, and, while it took many years to detect this feature 
(Cowie & Hu 1998; Hu et al. 1998), Lya emitting galaxies 
(LAEs) are now routinely observable from z ~ 0.2 (Deharveng 
et al. 2008; Cowie et al. 2010) to z > 7 (Hu et al. 2010; Ouchi 
et al. 2010; Lidman et al. 2012; Ono et al. 2012). However, 
while the detection of Lya in the high-redshift universe is 
relatively common, the physics of this emission is still not well 
understood. Since Lya is a resonance transition, it is likely that 
each photon scatters many times off intervening neutral material 
before escaping into intergalactic space. As a result, even a small 
amount of dust should extinguish the line, and indeed, only 
~25% of Lyman-break galaxies (LBGs) at z ~ 2-3 have enough 
Lya in emission to be classified as an LAE (Shapley et al. 2003). 
While it is possible for dusty galaxies to create an escape path 
for Lya via supernova-blown bubbles and/or exotic geometry 
(e.g., Verhamme et al. 2012) most analyses suggest that the 
LAE population as a whole is made up of young, low-mass, 
low-metallicity systems, possessing relatively little interstellar 
dust (e.g., Gawiser et al. 2007; Guaita et al. 2011). 

To date, most Lya emitters have been detected via deep 
narrow-band imaging with 4 m and 8 m class telescopes (e.g., 
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Gronwall et al. 2007 ; Ouchi et al. 2008). These surveys generally 
extend to low Lya luminosities and sample a wide range of the 
high-redshift galaxy luminosity function. Unfortunately, in the 
continuum, LAEs are usually quite faint, which makes studying 
their spectral energy distributions (SEDs) difficult. As a result, 
most of our knowledge about those physical properties which 
are encoded in the objects’ SEDs — information such as stellar 
mass, extinction, and population age — has come from stacking 
techniques (e.g., Gawiser et al. 2007; Guaita et al. 2011). These 
analyses only yield estimates for a “typical” LAE and may be 
subject to serious systematic biases associated with the stacking 
techniques (Vargas et al. 2014). Moreover, those few programs 
that have sought to measure the SEDs of individual LAEs (e.g., 
Finkelstein et al. 2009; Nilsson et al. 2011; Yuma et al. 2010; 
Nakajima et al. 2012; McLinden et al. 2014) have generally been 
restricted to very small numbers of objects. These efforts have 
been able to provide hints as to the range of properties exhibited 
by the population, but have been unable to probe the statistics 
of the entire LAE population. Thus, while we have some idea 
about the mass and dust content of “representative” LAEs, the 
distribution of physical parameters for the entire population 
remains poorly constrained. 

Here, we investigate the stellar populations of luminous 
Lya emitters by analyzing the individual SEDs of 63 1.9 < 
z < 3.6 LAEs detected by the McDonald 2.7 m telescope’s 
Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) 
Pilot Survey. In Section 2, we summarize the HETDEX Pilot 
Survey (HPS) and describe the ancillary ground-based, Hubble 
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Space Telescope ( HST) and Spitzer photometry which is avail- 
able for analysis. In Section 3, we briefly describe the SED- 
fitting code GalMC (Acquaviva et al. 2011) and the underlying 
assumptions used to derive stellar mass, extinction, and age 
from a set of broadband photometry which extends from the 
rest-frame UV through to the near-IR. We also outline the pro- 
cedures used to measure the physical sizes of the LAEs in a 
manner that is insensitive to the effects of cosmological surface 
brightness dimming. In Section 4, we present our results and 
show that the population of luminous z ~ 3 LAEs is quite het- 
erogeneous, with sizes extending from 0.5 kpc < r < 4 kpc, 
stellar masses ranging from 7.5 < logM/M Q < 10.5, and dif- 
ferential extinctions varying between 0.0 < E(B — V) < 0.4. 
We illustrate several trends involving LAE physical parameters, 
including a positive correlation between reddening and stellar 
mass, a positive correlation between stellar mass and galactic 
age, and a positive correlation between galaxy size and hya 
luminosity. We also examine the possible evolution of physical 
properties with redshift and compare our LAEs to other z ~ 3 
objects on the star-forming galaxy main sequence. We conclude 
by discussing the implications of our results for the underlying 
physical mechanisms of hya escape in high redshift galaxies. 

For this paper we adopt a cosmology with Hq = 
70 km s -1 Mpc -1 , Qm = 0.3, and = 0.7 (Planck 
Collaboration 2013; Hinshaw et al. 2013). 

2. OUR SAMPLE 

The LAEs chosen for study were discovered with the 
George & Cynthia Mitchell Spectrograph (previously known as 
VIRUS-P; Hill et al. 2008) on the 2.7 m Harlan J. Smith Tele- 
scope during the HPS (Adams et al. 2011). This integral-field 
instrument, which employs an array of 246 4" 2 diameter fibers, 
covers ~3 arcmin 2 of sky at a time, and delivers 5 A resolution 
spectra between the wavelengths 3500 A and 5800 A. The HPS 
itself surveyed a total of 169 arcmin 2 in the COSMOS (Scoville 
et al. 2007), GOODS-N (Giavalisco et al. 2004), MUNICS-S2 
(Drory et al. 2001), and XMM-LSS (Pierre et al. 2004) fields 
and reached a limiting line flux of 6.7 x 10 -17 erg cm -2 s -1 at 
5000 A (for 50% of its pointings) and 1.0 x 10 -16 erg cm -2 s -1 
at 5000 A (for 90% of the pointings). The final HPS catalog con- 
sists of coordinates, redshifts, R-band magnitudes, line fluxes, 
and equivalent widths for 397 emission-line selected galaxies. 
Ninty-nine of these sources are non X-ray emitting LAEs with 
1.9 < z < 3.8, rest- frame equivalent widths EWo > 20 A, and 
hya luminosities between 2.6 x 10 42 and 1.1 x 10 44 erg s 1 
(Adams et al. 2011; Blanc et al. 2011). A total of 74 of these 
HPS LAEs lie in the GOODS-N and COSMOS fields, where 
deep HST data is available. The redshift range for this subsample 
is 1.9 < z < 3.6. 

The process of assigning an optical counterpart to each HPS 
emission-line detection was challenging. As pointed out by 
Adams et al. (2011), there is an order of magnitude mismatch 
between the spatial resolution obtained from the 4" 2 diameter 
fibers of VIRUS-P, and that delivered by the broadband imagers 
of HST Thus, each assignment was done in a Bayesian manner, 
by calculating the likelihood of association for each object 
within a 10" window of the nominal position obtained from 
the spectroscopy (see Section 5.3 of Adams et al. 2011). 
Formally, the median probability for identifying the correct 
optical counterpart was 64%. However, as discussed below in the 
final paragraph of Section 3.1, there is no statistical difference 
between the distribution of SED properties for a sample LAEs 


with high-probability and/or confirmed counterparts and that 
for the sample of lower-probability associations. 

In total, we identified 67 HPS LAEs with optical counterparts. 
Four of these objects (HPS IDs 144, 145, 160, and 196) were 
removed from our analysis based on the work of Blanc et al. 
(2011), who showed that their UV slopes were more consistent 
with those of foreground [O ii] emitters than LAEs. This culling 
left us with 63 objects for analysis. Since the X-ray data in 
GOODS-N and COSMOS is deep enough to rule out most active 
galactic nuclei, we believe that the bulk of these objects are true 
hya emitting sources with 1.9 < z < 3.6 and monochromatic 
hya luminosities between 3.4 x 10 42 and 3.8 x 10 43 erg s 1 . 

Before proceeding further, we should note that the LAEs 
discovered by the HPS are significantly more luminous than 
the hya emitters found by most narrow-band surveys. While 
the 2 < z < 3 observations of Gronwall et al. (2007), Ouchi 
et al. (2008), and Guaita et al. (2010) typically reach hya 90% 
completeness levels of L(hya) ~ 10 42 erg s -1 , the median HPS 
limit is five times brighter than this. On the other hand, since 
the hya luminosity limit of the HPS is very nearly constant 
across the survey’s entire spectral range (see Figure 1 of Blanc 
et al. 2011), the data set covers an order of magnitude more 
co-moving volume than a typical narrow-band survey, with 
V = 5.63 x 10 5 Mpc 3 in the COSMOS and GOODS-N regions 
alone. This allows us to obtain good statistics on the bright end 
of the LAE population, and explore evolution over ~1.6 Gyr of 
cosmic time. 

The HPS fields, and in particular, the COSMOS and 
GOODS-N regions are rich in deep, ancillary imaging and pro- 
vide up to 18 photometric data points for SED fitting. Tables 1 
and 2 summarize these data. Most of the HPS/COSMOS and 
HPS/GOODS-N fields are part of CANDELS (Grogin et al. 
2011; Koekemoer et al. 2011), and 49 out of our 63 LAEs 
have deep HST optical and near-IR photometry (Song et al. 
2014; Finkelstein et al. 2013) from this program. Moreover, 
all of our targets have photometry from Adams et al. (2011), 
who used the deep ground-based images of COSMOS as their 
source frames (Capak et al. 2007). Note that many of the LAEs 
targeted in this survey are too faint in the continuum to be 
present in the published COSMOS and GOODS photometric 
catalogs; for these objects, AB magnitudes were determined 
by re-measuring the original images using the positions of the 
HST optical counterparts. Still, there are some non-detections. 
When this occurred, an upper flux limit was assigned as the 
lcr uncertainty of the local sky value. In some cases, these 
limits were crucial for constraining the SED properties of our 
targets. 

Data at longer wavelengths come from observations with the 
Spitzer telescope. Once again, most LAEs are far too faint to be 
present in the S-COSMOS and GOODS-N Spitzer catalogs, 
as these analyses have relatively high detection thresholds 
(1 /xJy in IRAC channel 1). Since the rest-frame near-IR is 
extremely important for determining stellar mass, we performed 
our own aperture photometry on these frames using MOPEX 8 
(Makovoz & Marleau 2005) at the known LAE positions. 
After experimenting with a variety of apertures, we settled 
on a photometric radius of 3"6, and then applied an aperture 
correction as described in the IRAC Instrument Handbook 9 . 


8 Information on MOPEX is available at 
http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/ 
tools/mopex/. 

9 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ 
iracinstrumenthandbook/28/ 
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Table 1 

COSMOS Field Photometry 


Hagen et al. 


Telescope 

Instrument 

Filter 

Central X 
(A) 

Original 

Survey 

Photometry 

5d Limits 
(AB) 

CFHT 

Megaprime 

w* 

4065 

COSMOS 

Adams et al. (2011) 

26.5 

Subaru 

Suprime-Cam 

B 

4788 

COSMOS 

Adams et al. (2011) 

27.4 

Subaru 

Suprime-Cam 

V 

5730 

COSMOS 

Adams et al. (2011) 

27.2 

Subaru 

Suprime-Cam 

r+ 

6600 

COSMOS 

Adams et al. (2011) 

26.9 

HST 

ACS 

F814W 

7461 

CANDELS 

Song et al. (2014) 

27.5 

Subaru 

Suprime-Cam 

/+ 

7850 

COSMOS 

Adams et al. (2011) 

26.9 

Subaru 

Suprime-Cam 

z+ 

8700 

COSMOS 

Adams et al. (2011) 

25.6 

UKIRT 

WFCAM 

J 

12850 

COSMOS 

Adams et al. (2011) 

23.6 

HST 

WFC3 

F125W 

13250 

CANDELS 

Song et al. (2014) 

26.4 

HST 

WFC3 

F160W 

14460 

CANDELS 

Song et al. (2014) 

26.5 

CFHT 

WIRCAM 

K 

21400 

COSMOS 

Adams et al. (2011) 

23.6 

Spitzer 

IRAC 

Channel 1 

37440 

S-COSMOS 

This paper 

23.9 

Spitzer 

IRAC 

Channel 2 

44510 

S-COSMOS 

This paper 

23.3 

Spitzer 

IRAC 

Channel 3 

59950 

S-COSMOS 

This paper 

21.3 

Spitzer 

IRAC 

Channel 4 

84870 

S-COSMOS 

This paper 

21.0 

Note. CANDELS covers 32 of 42 objects in this field. 







Table 2 





GOODS-N Field Photometry 



Telescope 

Instrument 

Filter 

Central X 

Original 

Photometry 

5(7 Limits 




(A) 

Survey 


(AB) 

Mayall 

MOSAIC 

U 

4065 

GOODS 

Adams et al. (2011) 

27.1 

HST 

ACS 

F435W 

4570 

CANDELS 

Finkelstein et al. (2013) 

27.8 

Subaru 

Suprime-Cam 

B 

4788 

GOODS 

Adams et al. (2011) 

26.9 

Subaru 

Suprime-Cam 

V 

5730 

GOODS 

Adams et al. (2011) 

26.8 

Subaru 

Suprime-Cam 

r+ 

6600 

GOODS 

Adams et al. (2011) 

26.6 

HST 

ACS 

F606W 

6690 

CANDELS 

Finkelstein et al. (2013) 

27.6 

HST 

ACS 

F775W 

7380 

CANDELS 

Finkelstein et al. (2013) 

27.5 

Subaru 

Suprime-Cam 

i+ 

7850 

GOODS 

Adams et al. (2011) 

25.6 

HST 

ACS 

F850LP 

8610 

CANDELS 

Finkelstein et al. (2013) 

27.3 

Subaru 

Suprime-Cam 

z+ 

8700 

GOODS 

Adams et al. (2011) 

25.4 

HST 

WFC3 

F105W 

11783 

CANDELS 

Finkelstein et al. (2013) 

26.6 

HST 

WFC3 

F125W 

13250 

CANDELS 

Finkelstein et al. (2013) 

26.4 

HST 

WFC3 

F160W 

14460 

CANDELS 

Finkelstein et al. (2013) 

26.5 

UH 2.2 m 

QUIRC 

H + K' 

20200 

GOODS 

Adams et al. (2011) 

22.1 

Spitzer 

IRAC 

Channel 1 

37440 

GOODS 

This paper 

23.9 

Spitzer 

IRAC 

Channel 2 

44510 

GOODS 

This paper 

23.3 

Spitzer 

IRAC 

Channel 3 

59950 

GOODS 

This paper 

21.3 

Spitzer 

IRAC 

Channel 4 

84870 

GOODS 

This paper 

21.0 


Note. CANDELS covers 17 of 21 objects in this field. 


3. DATA ANALYSIS 
3.1. Spectral Energy Distribution Fitting 

The SED of a galaxy encodes a number of physical param- 
eters, including stellar mass, age, dust content, and the current 
star formation rate (SFR). For example, since a galaxy’s near- 
IR flux arises principally from the evolved stars of all stellar 
populations, that part of the SED traces the system’s total stellar 
mass (Bell & de Jong 2001; Zibetti et al. 2009). In contrast, 
the slope of a galaxy’s far UV (~ 1600 A) continuum is fixed 
by the Rayleigh- Jeans tail of the blackbody emission from hot, 
young stars. The amplitude of the UV continuum thus yields 
the SFR and any flattening of the UV continuum’s slope is most 
likely due to the effects of dust (Kennicutt 1998; Calzetti 2001). 
Estimates of population age come primarily from the regions 
in between, as features such as the B aimer and 4000 A breaks 
are sensitive to the main sequence turnoff and the exact mix of 
intermediate age stars (e.g., Kauffmann et al. 2003). 


To extract this information, we began with the population 
synthesis models of Bruzual & Chariot (2003), which were 
updated in 2007 (BC07) with an improved treatment of the 
thermal-pulsing asymptotic giant branch (TP-AGB) phase of 
stellar evolution. This phase of stellar evolution can be impor- 
tant for systems older than ~10 8 yr, which is ~30% of our 
sample (see Section 4.4). We also performed fits using the older 
BC03 models, but due to the generally young ages of the stellar 
populations, these fits were statistically indistinguishable from 
the 2007 models. For the remainder of this paper, we will only 
refer to our BC07 results. For consistency with the works of 
Guaita et al. (2011) and Acquaviva et al. (2011), we adopted 
a Salpeter (1955) initial mass function (IMF) over the range 
0.1M o < M < 100 Mq, a Calzetti (2001) extinction law, 
and a Madau (1995) model for the effects of intervening inter- 
galactic absorption. Since stellar metal abundances are poorly 
constrained by broadband SED measurements, we fixed the 
metallicity of our models to Z = 0.2 Z Q ; this is roughly the 
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Observed Wavelength (//m) 



Figure 1 . Results of our SED fit to the photometry of HPS189. The blue points 
are the observed flux densities, the black line is the best fit SED, and the red 
triangles show the predicted flux density within each band. The left axis defines 
the scaling for the galactic continuum; the right axis gives the monochromatic 
flux scale applicable to the emission lines. 

(A color version of this figure is available in the online journal.) 

gas-phase abundance inferred from recent near-IR spectroscopy 
of z ~ 2 LAEs (Nakajima et al. 2013; Song et al. 2014). Emis- 
sion lines and nebular continuum, which can be an important 
contributor to the broadband SED of high-z galaxies (e.g., Atek 
et al. 2011; Schaerer & de Barros 2009), were modeled follow- 
ing the prescription of Acquaviva et al. (2011) with updated 
templates from V. Acquaviva (2012, private communication). 
Finally, following Guaita et al. (2011), we adopted the simple 
assumption that the SFRs of our LAEs have been constant with 
time. 

In keeping with these assumptions, we did not use any of the 
bandpasses listed in Tables 1 and 2 that lie redward of 3.3 /xm 
in the rest-frame. This is where the first polycyclic aromatic 
hydrocarbon line is located, and such interstellar medium (ISM) 
features are not accounted for in the stellar populations models. 
We also did not use data blueward of Lya, as the Madau (1995) 
correction is statistical in nature, and large excursions from the 
norm could bias our reddening measurements. 

Since SED fitting is a notoriously non-linear problem that 
may involve many local minima, highly non-Gaussian errors, 
and degeneracies between parameters, we chose to analyze 
our data using GalMC, a Markov-Chain Monte-Carlo code 
with a Metropolis-Hastings sampler (Acquaviva et al. 2011). 
This approach is much more computationally efficient than 
traditional grid searches, as it explores all regions of parameter 
space, while still spending the bulk of its time in the highest- 
likelihood parts of the probability distribution. The uncertainties 
associated with the fitted parameters are also much more 
realistic than those estimated using a simple x 2 minimization, 
as degeneracies between the variables are bettered explored and 
quantified. For each SED, three free parameters, stellar mass, 
E(B — V ), and age (which, under the assumption of a constant 
SFR history, is equivalent to SFR), were fit using four chains 
initiated from random starting locations. Once completed, the 
chains were analyzed via the CosmoMC program GetDi st (Lewis 
& Bridle 2002), and, since multiple chains were computed for 
each object, the Gelman & Rubin (1992) R statistic was used to 
test for convergence using the criterion R — 1 < 0. 1 (Brooks & 
Gelman 1998). 

As discussed at length by Conroy (2013), the results of 
our fits should be robust within the context of our model 
assumptions. Of course, any change to these assumptions will 


result in a systematic error. For example, the use of a Chabrier 
(2003) or Kroupa (2001) IMF would systematically reduce our 
stellar mass estimates by ~0.3 dex, while leaving our values 
for extinction and age essentially unchanged (Papovich et al. 
2011). Similarly, a different treatment of the TP-AGB phase may 
change the stellar mass estimates by up to ~0.3 dex (Zibetti et al. 
2009), while a shift to solar metallicity will generally increase 
our masses by ~0.1 dex. A full discussion of these systematic 
uncertainties is given by Conroy (2013). 

Figure 1 demonstrates the effectiveness of our fits by compar- 
ing the broadband photometry to the best fit SED for HPS 189. 
Since our data contain many bands of photometry, the SEDs of 
our galaxies are generally well constrained. Table 3 summarizes 
the SED-based properties of all 63 LAEs in our sample. 

As stated in Section 2, the assignment of optical counterparts 
to the HPS -discovered emission lines is probabilistic in nature. 
To investigate this further, we constructed a “clean” sample 
of HPS -LAEs, using a set of 29 candidates with either spec- 
troscopic confirmations (Finkelstein et al. 2011; Chonis et al. 
2013; Song et al. 2014) or a very high (^0.9) probability of 
association. We then used the Kolmogorov-Smirnov statistic to 
test whether the sample’s distributions of stellar mass, age, and 
reddening were in any way different from those formed from 
the remaining 34 LAEs. Figure 2 compares the empirical cu- 
mulative distribution functions (ECDFs) of the samples for all 
three SED parameters. In all cases, the two distributions are sta- 
tistically indistinguishable. Misidentifications are therefore not 
biasing the results of our SED fitting. 

3.2. Size Measurements 

To complement our SED-based estimates of stellar mass, 
extinction, and age, we also measured the sizes of the LAEs’ 
star-forming regions, via rest- frame UV measurements from 
HST. In the COSMOS field, deep F814W images are available, 
which, for the LAEs targeted in this study, sample the rest- frame 
wavelength region between ~2100 A and ~2800 A, depending 
on redshift. For GOODS-N LAEs, we have access to F606W 
observations, which are sensitive to rest- frame wavelengths 
from ~1600A to ~2100A. In both cases, we followed the 
exact same procedures described in Section 3.2 of Bond et al. 
(2009). After creating cutouts around each galaxy, we performed 
object identifications and background subtraction using the 
routines found in SExtractor (Bertin & Arnouts 1996). We then 
obtained a measure of size by using the IRAF 10 program phot 
to determine each object’s flux- weight centroid and magnitude 
through a series of circular apertures. These aperture magnitudes 
were then used to define each LAE’s half-light radius. As 
described by Bond et al. (2012), the uncertainty in this type 
of measurement is given by 

— = 0.54 —t, (1) 

r f 

where r is the half-light radius, / and erf are the flux and 
associated flux uncertainty, and o r is the resultant error on 
the half-light radius. For our sample of LAEs as a whole, the 
fractional median uncertainty for the measured half-light radius 
is 4%. We note that this measure of size is much less sensitive to 
the (1 + z) 4 effects of cosmological surface brightness dimming 
than estimates based on limiting isophotes. Furthermore, the 

10 IRAF is distributed by the National Optical Astronomy Observatory, which 
are operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
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Table 3 

SED Fitting Results 


Galaxy ID 

Redshift 

Log Stellar Mass 

Log Age 

£(£ - V) 

Half-light Radius 

Xv 



(M 0 ) 

(yr) 


(kpc) 


HPS150 

2.90 

Q oc+0.14 

7 cq+0.30 
'• Jy -0.2S 

Oil 8+0015 
U - A 1 °- 0.012 

1.19 ±0.02 

4.1 

HPS153 

2.71 

9 50 +0 - 20 
y,JU -0.19 

7-43 !° 0 5 4 2 8 

0 289 +0 - 026 

O.2O7_q 016 

1.62 ±0.06 

1.5 

HPS154 

2.87 

Q 07+0.34 
7 -0.23 

6.87+ 1 !™ 

0 747 + 0.059 
U.7HO_oo07 

0.52 ±0.11 

3.5 

HPS161 

3.25 


7 70 + 0.23 
7 - Z °-0.24 

0 775 + 0.012 
kj.j /J _ 0.007 

2.01 ±0.06 

1.0 

HPS164 

2.45 

Q Q7+0.13 
y ’ y/ -0.14 

7 87+0.28 
7 -° 7 -0.30 

0 7O4 +0 - 019 
U.7U4 -_ 0 012 

1.96 ±0.05 

2.9 

HPS168 

3.45 

0 44 +O .25 
5 - 44 -0.24 

6 Q7+0.65 
7 -0.74 

0 O7O +0 020 
U.U /u _ 0 070 

1.53 ±0.18 

2.8 

HPS174 

3.45 

0 77 +O .59 
^• 72 -0.51 

7 67+ 1-00 

7 7 -0.98 

H 778+0-037 
KJ ’ ZjD °-0.238 

0.87 ±0.28 

3.2 

HPS182 

2.43 

9.18t° 0 3 3 4 6 

7 88+0-84 
7 '°°-0.80 

0 147 + 0.048 

U. l^7_o 443 

0.82 ±0.09 

5.0 

HPS183 

2.16 

8-64!°o 3 3 0 2 

7-8® 

0 052 +0 - 015 
u.ujz _ 0 Q52 

1.25 ±0.10 

2.3 

HPS184 

3.21 

Q oq+0.27 
.27 

8 - 20 t ° 0 7 7 5 2 

0 17O +0 - 050 
U * 1 /u -0.170 

0.72 ±0.14 

1.4 

HPS189 

2.45 

8.64-° d 1 1 8 9 

7 63 +0 - 33 
7 — 0.34 

0 O72 +0 - 014 

U.U 7 -7_ 0012 

0.94 ± 0.06 

1.0 

HPS194 

2.29 

9 28 +0 - 28 
^• Z °-0.33 

8 1 /:+0.47 
-0.59 

0 O46 +0 015 
u.u^o _ 0046 

1.96 ±0.07 

7.9 

HPS197 

2.44 

7-83^, 

6 O 1 + 0.88 
°* yA - 1.01 

0-093*™ 22 

0.57 ±0.11 

1.1 

HPS205 

2.91 

8-™t°o 4 5 8 o 

8 77 + 0.66 
^ 7 - 0.68 

0 . 068 !°(j 26 2 8 

1.55 ±0.53 

6.2 

HPS207 

2.71 

7 qq+0.36 
7 -^-0.38 

7 aq+0.64 
7 •■ 2 — 0.60 

0 045 +0 - 01 ° 
u.u+j _ 0Q45 

0.97 ±0.17 

2.8 

HPS210 

3.49 

0 70 +O. 3 I 
°- Z3 -0.31 

7 26 +0 - 50 

7 - zo _o.47 

0 01 7 +0-003 
U.U17 — 0 013 

1.59 ±0.26 

3.1 

HPS213 

3.30 

Q QO+018 
y * yU - 0.21 

9 01 +0 - 27 
y ' U1 -0.25 

0 O2O +0 004 

u.uzu _ 0 020 

1.39 ±0.06 

9.3 

HPS214 

3.30 

7 qc+0.53 
'• yj -0 .44 

6 04 +I.O 9 
°*^- 1.12 

0 125 +0 - 028 

U>1Z -0.125 

0.79 ± 0.25 

0.4 

HPS223 

2.31 

0 1 q+ 0.36 
°- iy -0.36 

7 Q4+0.50 
0.49 

0 029 + ° 004 
u.uzy _ 0029 

3.61 ±0.54 

3.2 

HPS229 

3.04 

Q 8 O +011 

0 41 +O.I 5 
o *^-0.14 

0 01 5 + 000 3 

U.Ul J _ 0 015 

1.33 ±0.03 

6.3 

HPS231 

2.72 

7 RQ + 0 - 27 
/ - oy -0.28 

V-657°d1g, 

0 03 1 + 0 - 006 

U.U 7 i _ 0 031 

1.63 ±0.21 

2.1 

HPS244 

2.10 

7 64 +0 - 15 

6-38t° 0 5 7 7 4 

0 O7O +0 014 
u.u/u _ 0070 

1.05 ±0.16 

0.5 

HPS249 

3.27 

9.06^ 

7 62 +0,38 
7 • oz — 0.37 

0 266 +0 - 019 
u.zoo _ 0018 

1.43 ±0.17 

3.0 

HPS251 

2.29 

8-4® 

7 C7+0.40 

7 ^ 7 -0.42 

0 025 +0 - 006 
u.uzd _ 0025 

1.04 ±0.05 

4.4 

HPS253 

3.18 

001 +0.25 
-0.25 

7 60 +0 - 35 

7 -0.35 

0 020 + ° 005 
u.uzu _ 0020 

2.07 ±0.10 

2.1 

HPS256 

2.49 

0 17 +O.I 9 
7 — 0.18 

7 70 +O. 3 I 
7 - Z °-0.32 

0 041 +0 01 ° 

U.U41-0 041 

1.07 ±0.07 

2.9 

HPS258 

2.81 

0 04 +O.I 7 
6 ‘ y4 -0.16 

7 87+0-26 
7 -0.25 

0 022 +0 - 005 
u.uzz _ 0 022 

2.57 ± 0.06 

7.3 

HPS263 

2.43 

0 00 +O.I 8 

5.58_0 17 

7 60 +0 - 34 
7 — 0.31 

0 049 +0 012 

U.U^y_0049 

1.51 ±0.05 

3.8 

HPS269 

2.57 

0 07 +O.I 7 
°-° 7 -0.17 

0 7C+0.27 
0,z — 0.30 

0 025 +0 - 008 
u.ujj_o 035 

0.90 ± 0.07 

3.0 

HPS274 

2.87 

9 06 +012 
y * UO - 0.12 

7 45 + 0.21 
7 - 0.20 

0 120 + ° 009 
u.izu _ 0007 

1.01 ± 0.02 

4.3 

HPS283 

3.30 

g /:c+0.16 

8 59+0.24 
°* jy - 0.21 

0 027 +0 - 006 
u.uz / _ Q Q27 

2.31 ±0.10 

2.8 

HPS286 

2.23 

0 04 +O .25 
5 '^-0.28 

8 11 +0.37 
11 -0.41 

0 024 +0 - 004 
u.uz ^_ 0024 

1.97 ±0.10 

24.3 

HPS287 

3.32 

0 77 +O. 3 O 
7 7 -0.32 

7 47+0.48 
7 7 -0.48 

0 2O7 +0 - 020 

U.Z.U / —o 019 

0.37 ±0.21 

1.3 

HPS288 

3.04 

0 77 +O .28 
7 7 -0.29 

7 71 +0.48 
/ '- 31 -0.37 

0 O29 +0 01 ° 
U.U7y_ Q 039 

1.35 ±0.04 

0.8 

HPS292 

2.87 

8-59^0*20 

771 +0.28 
7 - 71 -0.29 

0 022 +0 - 006 
u.uzz_o 022 

0.73 ± 0.04 

5.8 

HPS296 

2.84 

8 - 2 ® 

7 ci +0.30 
7 *^ — 0.31 

0 01 i+o o° 2 
U - Uii -0.011 

0.97 ±0.11 

9.4 

HPS306 

2.44 

9 01 +012 
y ' Ui - 0.12 

7 77+0.21 
7 -0.20 

0 079 +0 009 

U.U77_q 039 

1.35 ±0.04 

3.3 

HPS310 

3.07 

9 26 +0 - 22 
y - ZO -0.18 

0 C 4 +O .32 
^-0.25 

0 024 +0 - 004 
U.UZ ‘+_ 0 024 

0.87 ±0.05 

4.1 

HPS313 

2.10 

Q O 7 +O.O 8 
^• oz — 0.08 

7 74 +O.H 
7 * 74 - 0.12 

0 161 +0007 
U - i0i -0.006 

2.34 ± 0.02 

17.3 

HPS315 

3.07 

Q O4 +0 - 23 
y ' U4 -0.23 

7 49+0.37 
7 * 4 Z -0.34 

0 O96 +0 011 

u.uyo_ 0 010 

3.03 ± 0.09 

2.0 

HPS316 

2.81 

o 40 +O.H 
y '^ u - 0.10 

0 78+0.19 
°- Z °-0.18 

0 027 +0 - 006 
u.uz / _ 0 027 

1.26 ±0.03 

1.9 

HPS318 

2.46 

Q 54+0-12 

7 7Q+0.23 
7 ' 7 —0.23 

0 1 79 +O.OH 

U. 17^-0 010 

2.42 ± 0.07 

6.9 

HPS338 

2.60 

7 86+0.43 
'•° b -0.30 

6 - 8 ® 

0 1 1 6 +0014 
u * iiO -0.116 

3.52 ±0.13 

0.9 

HPS341 

2.93 

0 41 +0.18 
5 - 41 -0.19 

6-91!°o 5 5 5 5 

0 059 +0 - 017 
u.ujy _ 0059 

2.53 ±0.50 

1.0 

HPS360 

2.92 

9-8® 

7 96+0.56 
7 — 0.53 

0 266 +0 046 
u.zoo _ 0017 

1.61 ±0.04 

1.9 

HPS370 

3.18 

8 . 40^1 

7 78+0-33 
L ^- 0.32 

0 022 +0 - 008 
U.Ujj_q 033 

0.94 ± 0.02 

0.9 

HPS372 

2.76 

7 C 7 + 0.44 
7 *^ Z — 0.39 

7 O 8 +0 ' 97 
7 * u ^_ 1 18 

0 O67 +0 - 017 
U.UO /-0 067 

0.36 ± 0.06 

2.9 

HPS389 

2.59 

9 09 +01 ° 
y - uy -o.io 

7 77 +O.I 7 
7 * /Z -0.17 

0 107 +001 ° 
U - iU/ -0.008 

1.28 ± 0.02 

1.9 

HPS391 

2.96 

Q 44+0.13 
*•^-0.15 

8-132°6 2 3 7 3 

0 ioi +0 016 

U>1U1 -0.015 

0.96 ±0.01 

1.2 

HPS395 

2.27 

Q 74+0.09 
y,:) ^-0.08 

8 . 36 ^o!i 7 

0 O26 +0 009 
U.U70_o 036 

3.79 ±0.52 

1.6 

HPS402 

2.97 

8.36^042 

6 97 + 0.28 
7 -0.21 

0 05 1 + 0 - 008 

u.uj i_o 008 

1.87 ±0.03 

2.2 

HPS403 

3.18 

8 - 9 ®> 

0 1 7+0.41 
°- A 7 -0.46 

0 027 +0 - 010 
U.UJ 7 -0.037 

2.70 ±0.16 

7.8 

HPS415 

3.37 

q 71 + 0.15 
y - 71 -0.18 

9 05 + 0.22 
y,U - 0.21 

0 01 9+000 4 

U.Uiy_ooi9 

1.09 ±0.03 

4.3 
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Table 3 

(Continued) 


Galaxy ID 

Redshift 

Log Stellar Mass 
( Mq ) 

Log Age 
(yr) 

E(B - V ) 

Half-light Radius 
(kpc) 

Xv 2 

HPS419 

2.23 

1 O.OO+q ‘jj 

0 oo+0.28 
5.53—0 25 

0 1 82 +0 - 019 
0.014 

1.67 

± 

0.02 

1.8 

HPS420 

2.93 

8-6V+°o 1 1 8 8 

771 +0.28 
' ’ ' —0.28 

0 052 +0 - 011 
O.OJZ— 0052 

0.56 

± 

0.01 

5.4 

HPS426 

3.40 

Q C7+0.20 
y ' J / - 0.20 


0 049 +0 - 011 

U.U^V— 0049 

0.70 

± 

0.02 

4.1 

HPS428 

3.34 

Q 46+0.10 
y - 4O -0.11 

8 -04-o.20 

0 O72 +0 - 010 
U.U /Z— 0 009 

1.93 

± 

0.07 

2.9 

HPS434 

2.27 

9.03^ 

0 74+0.71 
/4 -0.56 

0 078 +0 026 

U.U ' 3—0 073 

1.03 

± 

0.04 

3.1 

HPS436 

2.42 

8-36+S 

7 70+O.24 
' '^-0.26 

0 O77 +0 - 009 
U - U/ ' -0.007 

1.78 

± 

0.02 

1.5 

HPS447 

3.13 

8.39+Oj‘o 

6 76+ 0 - 60 
°- :>O -0.74 

0 022 +0 009 
U.U33— 0033 

3.60 

± 

0.11 

2.5 

HPS462 

2.21 

lO+S+O; 1 , 3 , 

0 f . 1+0.32 
5.01-0 29 

0 296 +0 - 024 
U.ZVO— 0 020 

1.38 

± 

0.84 

2.8 

HPS466 

3.24 

q 04+ 0 03 
y ' U4 -0.03 

z 07+O.O5 
°-° Z -0.04 

0 1 74+ 0 - 004 

U. 1^^—0 004 

1.51 

± 

0.01 

12.9 

HPS474 

2.27 

0 Q7+O.2O 
°- yi - 0.20 

7 Q5+0.35 
'• yj - 0.33 

0 096 +0 - 014 
U.UVO-0010 

1.18 

± 

0.02 

8.3 


observations differ in depth (COSMOS is a single orbit while 
GOODS is 2.5 orbits), and so any biases from surface brightness 
limits can be found by comparing the half-light radii from both 
surveys. A two-sample K-S test showed that the distributions 
of half-light radii derived from GOODS and COSMOS are not 
significantly different and thus we should not be concerned 
about effects from cosmological surface brightness dimming. 

4. RESULTS 

Table 3 gives the best-fit solutions to our SED fits, their 
reduced x 2 values, and the LAE’s half-light radii as measured 
on the HST frames. We discuss these results below. 

4.1. Size 

The left-hand panel of Figure 3 shows the distribution of half- 
light radii for the 63 luminous Lya emitters in our sample. 11 
From the figure, it is clear that the highly luminous LAEs of 
the HPS have a wide range of sizes: while the peak of the 
distribution is close to 1 .2 kpc, there is a distinct tail that extends 
all the way out to ~4 kpc. The median of the distribution is 
1.35^0% kpc, where the 68% confidence interval comes from a 
bootstrap analysis (Efron 1987). For comparison, the typical size 
of the narrow-band selected 2 < z < 3 LAEs studied by Bond 
et al. (2012) is ~1.0 kpc, while that for z ~ 3 LBGs is close 
to 4 kpc (Pentericci et al. 2010). As demonstrated by Spearman 
rank correlation coefficient (p = —0.15) and Figure 4, these 
sizes show no significant correlation with redshift; this result 
is consistent with that of Malhotra et al. (2012), who also saw 
no size evolution in samples of narrow-band selected LAEs 
between 2.5 < z < 6. This is in contrast with the strong 
evolution seen in the sizes of LBGs over the same redshift 
range (e.g., Ferguson et al. 2004). 

On the other hand, as the right-hand panel of Figure 3 shows, 
there is a relation between LAE half-light radius, as measured in 
the rest- frame UV, and the amount of luminosity emitted in the 
Ly a line. Although the scatter in the diagram is substantial, the 
Spearman rank order coefficient reveals a positive correlation 


1 1 This type of plot will be used throughout this work and shows a histogram, 
a kernel density estimation (KDE), and the empirical cumulative distribution 
function (ECDF). All KDEs contained herein use a Gaussian kernel and were 
calculated using the density function in R (R Core Team 2012). The bandwidth 
for each KDE was found using the rule of thumb from Scott (1992). A simple 
change in the choice of bins for a histogram can change the interpretation of 
the science; since KDEs do not require binning, they do not suffer from this 
effect. Every ECDF plotted will also have dotted upper and lower limits, which 
represent the 1 a asymptotic errors (Donsker 1952). 


between the two variables (p = +0.31), which is significant 
at the 2.5a (99%) level. Moreover, up to a half-light radius of 
2.5 kpc, p = +0.45, implying a 3.4a (99.9% confidence) result. 
Within this range, Ly a luminosity appears to grow linearly with 
galaxy size with a slope of 0.6 zb 0.2 x 10 43 erg s -1 kpc -1 . 
Unfortunately, the data outside this range are too sparsely 
populated to draw any conclusions. 

4.2. Stellar Mass 

At 2 < z < 3, SED stacking analyses have produced 
estimates for the median LAE mass in the range of ~10 8 to 
10 10 Mq, depending on whether the samples under study were 
detectable on Spitzer/ IRAC images (Lai et al. 2008; Guaita 
et al. 2011; Acquaviva et al. 2011). In the left-hand panel of 
Figure 5, we show the distribution of individual LAE masses 
derived from our constant SFR model. Although the median 
mass of \og(M / Mq) = 8.97+ lies in the range inferred 
from the previous stacking analyses, the data span a factor of 
a thousand, from M ~ 10 7 5 to 10 10,5 Mq. This distribution is 
consistent with that recently found by McLinden et al. (2014) 
for a set of extremely luminous (L(Lya) > 10 43 erg s l ) 
z ~ 3.1 LAEs found via narrow-band imaging. Interestingly 
if we fit our LAE mass distribution with a power law over the 
range 8.5 < lo gM/M© < 10.5, then the most-likely slope, 
a = — 1.3 =b 0.1, is similar, or perhaps only slightly shallower 
than the slope usually adopted for the epoch’s galaxy mass 
function (e.g., Muzzin et al. 2013; Tomczak et al. 2013). This 
suggests that bright Ly a emitters are drawn from an underlying 
distribution that is not strongly dependent on stellar mass. 
Moreover, as the right-hand panel of Figure 5 demonstrates, 
there is no obvious correlation between stellar mass and Lya 
luminosity: at any logL, one can find LAEs of all masses, 
and galaxies of any given mass can have a wide range of Ly ot 
luminosity. Unless this behavior abruptly changes at low Lya 
luminosity, it would appear that large- volume LAE surveys are 
an excellent way of sampling virtually the entire range of the 
high-redshift galaxy mass function. 

One additional feature of Figure 5 worth noting is the absence 
of LAEs with masses less than ~5 x 10 7 Mq. There are two 
possible reasons for this. The first is the depth of the imaging: 
there are seven sources in the HPS survey for which Adams 
et al. (2011) could find no obvious counterpart in the rest- frame 
UV. An examination of the CANDELS frames and grism spectra 
from the 3D-HST program (van Dokkum et al. 2013) reveals that 
only one of these missing objects has any detectable flux in the 
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log Age 



Figure 2. Empirical cumulative distribution functions of the clean sample of LAEs (blue), compared to those for LAEs with lower-probability counterparts (green). 
The dotted lines represent the 1 a asymptotic errors (Donsker 1952). In all cases the K-S test cannot reject the null hypothesis that both samples are drawn from the 
same underlying distribution. 

(A color version of this figure is available in the online journal.) 




Half-light Radius (kpc) 


Figure 3. Left: a histogram showing the distribution of half-light radii for the HPS LAEs. The cumulative distribution is shown in blue, while the kernel density 
estimation (KDE) of the distribution is in black. The dashed line illustrates the resolution limit of the data. The median half-light radius of 1.35 kpc is consistent 
with that found for narrow-band-selected LAEs at z ~ 2.1, but larger than the median at z ~ 3.1 (Bond et al. 2012). Right: a comparison of half-light radius to Lya 
luminosity. The gray points are the individual measurements, while the dark blue squares show the median value of each bin. The error bars in x delineate the size of 
each bin, while the error bars in y report the bin’s standard error. Although the scatter is large, the correlation for LAEs with r < 2.5 kpc is significant at the 3.4cr 
(99.9%) confidence level, and for the entire sample, the trend is confirmed with 2.5cr (99%) confidence. 

(A color version of this figure is available in the online journal.) 


rest-frame optical. (The counterpart of HPS 266 is detected at 
a (2000) =10:00:29.818, 5(2000) =+2:18:49.20.) The missing 
6 objects may therefore be part of the extreme low-mass tail of 
the mass function. 

A second possible explanation for the missing lower-mass 
galaxies comes from the limitations imposed by our input 
physics. Formally, the population synthesis models of Bruzual 
& Chariot (2003) are applicable to stellar systems with ages 
between 10 5 yr < t < 2 x 10 10 yr. However, the CB07 (and 
BC03) models used in our analysis, and the Padova isochrones 
upon which they are based, were developed using stars with ages 


of 10 6 ' 6 yr and older (Conroy 2013). Any system younger than 
this must therefore be subject to greater systematic uncertainties. 
The lack of systems with masses below x 10 7 M 0 may 
therefore be an artifact of our SED fitting. Furthermore, since 
our SED fits assume a constant SFR, very low mass systems 
almost certainly have very young ages, further increasing the 
potential for systematic errors at the very faint end of our sample. 

Figure 6 plots our derived stellar masses against two param- 
eters, LAE size and redshift. Neither shows a significant trend. 
While Bond et al. (2012) did find a correlation between mass and 
half-light radius, their analysis dealt with stacked images, not 
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Figure 4. Lack of LAE size evolution with redshift. The gray points are individual measurements while the blue squares are binned medians. A bias in size measurement 
from cosmic surface brightness dimming would manifest as a decrease in half-light radius as (1 + z) 4 . There is no evidence for this effect. 

(A color version of this figure is available in the online journal.) 



Figure 5. Left: the distribution of masses for the HPS LAEs under the assumption of a constant star formation rate throughout history. The items plotted on the left 
figure are described in the caption of Figure 3. The distribution is very nearly flat in the log over three orders of magnitude. Right: a comparison of galaxy mass with 
Lyc x luminosity. There is no significant correlation between the two parameters, implying that the sample’s Lya flux-limit does not propagate strongly into a constraint 
on stellar mass. 

(A color version of this figure is available in the online journal.) 


individual SEDs. Our null result also agrees with that found for 
surveys of UV-bright galaxies in the same redshift range. LAE 
sizes (as measured in the UV continuum) and stellar masses do 
not seem to be related (e.g., Mosleh et al. 201 1). 

4.3. E(B - V ) 

Since Lya photons resonantly scatter off interstellar hydro- 
gen, even a small amount of extinction can reduce the emergent 
emission-line flux by several orders of magnitude. Thus, it has 
long been argued that LAEs are mostly young, metal-poor ob- 
jects with very little dust (e.g., Gawiser et al. 2006; Mao et al. 
2007). Nevertheless, evidence for the existence of dust within 
LAEs has been seen in the work of Linkelstein et al. (2009) 
among others. As the left-hand panel of Ligure 7 illustrates, our 


data demonstrate that, indeed, most LAEs are dust-poor. Based 
on our SED analyses, half of the HPS LAEs have internal stel- 
lar reddenings E(B — V) < 0.07, though there is a tail that 
extends all the way out to E(B — V) ~ 0.4. The median of 
the E(B — V ) distribution is 0.067+ Q^g. Notably, all the high- 
extinction objects are drawn from the high-mass end of the LAE 
mass function: every LAE with E(B — V) > 0.25 has a mass 
greater than MO 9 M 0 . This agrees with correlations between 
mass and extinction seen in both the local universe and at high 
redshift (e.g., Garn & Best 2010; Kashino et al. 2013). 

Figure 8 displays our estimates of differential reddening 
versus redshift and galaxy size. There are no significant trends in 
either diagram. Unlike Guaita et al. (201 1) and Acquaviva et al. 
(20 1 1 ), we see no evidence for any change in the mean reddening 
of the LAE population versus redshift, and, unlike Bond et al. 
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Figure 6. Left panel displays stellar masses vs. redshift, with the individual LAEs shown in gray and the median values for each redshift bin plotted as blue squares. 
The right panel similarly compares stellar mass to half-light radius. The dashed line illustrates the limit of our spatial resolution. In both figures, the error bars in x 
illustrate the sizes of the bins, while the errors in y report the standard error within each bin. The black circles are averages of various continuum selected galaxies 
from Mosleh et al. (201 1) and the black triangles are results from stacks by Bond et al. (2012). 

(A color version of this figure is available in the online journal.) 



E(B-V) 



log Mass (M 0 ) 


Figure 7. Left: the distribution of stellar reddenings derived from our constant SFR models. The items plotted are described in the caption of Figure 3. The median of 
the distribution is E(B — V) = 0.07, but there is a tail that extends out to E{B — V) ~ 0.4. Right: the stellar reddening of our LAEs as a function of galaxy mass. 
The black circles represent the results of stacking analyses for LAEs at z ~ 2.1 and z ~ 3.1 (Acquaviva et al. 2011, 2012; Guaita et al. 2011). Low-mass objects are 
uniformly dust-poor, but objects with M > 10 9 M 0 can have a wide range of internal extinction. Stellar mass and extinction are correlated with a Spearman rank 
coefficient of p = 0.3, indicating 2.4cr (98.5%) significance. 

(A color version of this figure is available in the online journal.) 
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Figure 8. Comparisons of stellar reddening with redshift (left) and half-light radius (right). The gray points show individual LAEs while the blue points illustrate the 
median values within each bin. The dashed line shows the limit of our spatial resolution. Neither diagram displays a significant trend. 

(A color version of this figure is available in the online journal.) 
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Figure 9. Left: the age distribution of the HPS LAEs. The ages extend over ~2.5 dex. Right: the distribution of LAE ages vs. redshift, with individual objects plotted 
in gray and binned medians in blue. There is no evidence for evolution in the LAE population. 

(A color version of this figure is available in the online journal.) 



Figure 10. Dust- corrected UV star formation rate plotted against stellar mass 
for various high-redshift galaxies. Our luminous, spectroscopically selected 
LAEs are shown in red, and narrow-band-selected LAEs from Vargas et al. 
(2014) are in pink. The blue and green points show higher-mass galaxies from 
the Rodighiero et al. (2011) survey of COSMOS and GOODS-N: blue are 
continuum-selected ( BzK ) galaxies, while green shows submillimeter bright 
systems found by the Herschel - PACS instrument. The solid black line is the 
star-forming “main sequence” defined by Daddi et al. (2007) for 1.5 < z < 2.5; 
the extrapolation of this line to lower SFRs is shown with a dashed line. LAEs 
and submillimeter galaxies lie above this relation, along the main sequences 
found by Whitaker et al. (2012) for z = 2.0 (lower dotted line) and z = 3.5 
(upper dotted line). 

(A color version of this figure is available in the online journal.) 

(2012), we see no correlation between E(B — V) and half-light 
radius. Again, we caution that the LAE samples considered 
here are more luminous than those derived from narrow- 
band surveys, and by studying individual, rather than stacked 
spectra, we are avoiding many of the systematic difficulties that 
complicate the interpretation of previous measurements (Vargas 
et al. 2014). 

4.4. Age 

Figure 9 displays the age distribution for the HPS LAEs, 
under the simplest assumption of a constant SFR history. The 
two figures together support the stacking analyses of Lai et al. 


(2008), Guaita et al. (201 1), and Acquaviva et al. (201 1), which 
argued that LAEs are relatively young, with ages between 
~10 7 and 10 9 yr. The median age of the HPS sample is 
1 ogt = 7.96t°Q l i 4 , and just ~3% have ages greater than a Gyr. 
Interestingly, there no evidence for evolution in the sample. This 
disagrees with result of Guaita et al. (201 1), who found narrow- 
band selected z = 2.06 LAEs to be older and dustier than their 
z = 3.1 counterparts. It is also in conflict with the re-analysis by 
Acquaviva et al. (2012), who concluded that these same LAEs 
were “growing younger” with time. A likely explanation for this 
discrepancy lies in the details of the stacking procedure used by 
both groups, as slight differences can produce discrepant results 
(see Vargas et al. 2014). Our analysis of individual LAEs avoids 
that pitfall. 

4.5. Star Formation Rates and the Main Sequence of Galaxies 

The observed SFR of an LAE can be derived several ways: 
from the luminosity of its Ly a emission line via the assumption 
of Case B recombination, 

SFR(Lya) = 9. 1 x 10“ 43 L(Lya) M 0 yr -1 (2) 

(Brocklehurst 1971; Kennicutt 1998), from the extinction- 
corrected flux density of the UV continuum between 1500 A 
and 2800 A, 

SFR(UV) = 1.4 x 10“ 28 L y M q yr -1 (3) 

(Kennicutt 1998), and via the fit to the object’s SED (a value 
which is largely dependent on the UV emission, but which may 
also include factors such as age). 

Figure 10 plots our dust-corrected UV SFRs against stellar 
mass, and compares these data to those obtained from other 
samples of high-z galaxies. From the figure, it is clear that each 
selection technique identifies galaxies in a different region of 
the mass-SFR diagram. LAEs are primary low-mass, low-SFR 
objects that lie above the star-forming main sequence found by 
Daddi et al. (2007), in a region of the diagram consistent with 
measurements of galaxy main- sequence evolution (Whitaker 
et al. 2012). Both IFU and narrow-band- selected LAEs fall in 
this same region, confirming that both discovery techniques 
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Figure 11. Left: a histogram of the ^-parameters for our spectroscopically selected LAEs. Right: these same q values as a function of E(B — V). The data are biased, 
in that systems with large values of q will have Lya suppressed below the limit of detectability. All dusty systems have q near 1, while systems with small E(B — V) 
can have a range of q values. 

(A color version of this figure is available in the online journal.) 




Figure 12. Left: the escape fraction of Lya photons as a function of stellar mass. Note the negative correlation. Values greater than one can be due to an anisotropic 
ISM, as Ly a may be escaping through voids blown out by supernovae or other feedback. The decline in the maximum escape fraction with stellar mass supports 
this model, as it would take much more energy to blow out a hole in a massive galaxy. Right: the correlation between half-light radius and q\ galaxies with small or 
negative q factors are physically smaller. This is again consistent with the idea with anisotropic Lya emission, as it would be easier for supernovae to create holes in 
smaller, less massive galaxies. 

(A color version of this figure is available in the online journal.) 


trace the same population. UV continuum-detected ( BzK ) 
galaxies are higher-mass, high-SFR objects, while Herschel- 
PACS-selected submillimeter galaxies are high-mass objects 
that, like LAEs, fall predominantly above the star-forming main 
sequence (Rodighiero et al. 2011). As the distributions of LAEs 
and submillimeter galaxies abut each other, it is tempting to 
associate the two classes. If submillimeter systems are the results 
of merger-driven starbursts (Conselice et al. 2003), then LAEs 
could potentially be their low-mass and low-dust counterparts: 
in general, a starburst event will move a galaxy up and slightly 
to the right on this diagram. Gronwall et al. (2011), however, 
see no strong evidence for interactions in LAE morphologies. 

4.6. Lyot Escape Fractions and the q-f actor 

Like the UV flux density, the Ha emission line is a well- 
known and well-understood SFR indicator (Kennicutt 1998). 
Since under Case B recombination, roughly three quarters of all 
B aimer transitions produce a Lya photon, this means that Lya 


should also be a robust tracer of star formation. Thus, if all the 
Lya and UV continuum photons escape into intergalactic space, 
Lya and the dust-corrected UV continuum should be well- 
correlated. Systematic deviations from a one-to-one relation 
then measure the escape fraction of Lya photons. Note that 
this differential procedure sidesteps the issue of whether the 
reddening derived from the stars is the same as that for the 
gas, but it does assume that both the UV and Lya emission is 
isotropic. It is also susceptible to a systematic error associated 
with the timescale of star formation. The nominal conversion 
between UV emission and SFR (Equation (3)) assumes a 
timescale for star formation that is 10 times longer than that 
for emission line gas. If a system is undergoing a rapid burst of 
star formation (t < 10 8 yr), its two SFRs indicators may not 
be comparable. Nevertheless, the ratio of UV flux to Lya can 
provide constraints on the radiative transfer of the emission line. 

The question of the dependence of Lya escape fraction on 
SFR and redshift has been recently discussed in Dijkstra & 
Jeeson-Daniel (2013). This paper presents two models: one in 
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which the escape fraction is independent of SFR, and a second 
where the escape fraction decreases as the SFR increases. We 
find a significant inverse correlation between stellar mass and 
escape fraction, p = —0.54 (4.5a or >99.999% significance) 
which supports the second model; this is shown in Figure 12. 
Unsurprisingly, the Ly a escape fraction also correlates with 
differential extinction, as mass and E(B — V) are coupled (see 
Figure 7). We note that the median escape fraction of our sample, 
~0.5 (or 0.6, if we use the SED-based SFRs), is somewhat larger 
than the 0.29 value found by Blanc et al. (201 1) using the same 
sample of LAEs. Most of this difference is due to the use of 
deeper CANDELS data, which greatly improves the photometry 
and fixes the slope of the rest- frame UV. 

Perhaps a more useful way of looking at the radiative transfer 
problem is through the variable q , which relates Lya optical 
depth to that of the stellar continuum at 1216 A. As defined by 
Finkelstein et al. (2008), 


h — w 

U216 

where = 0.921 kxE(B — V ) and k\ 2 \e = 5.27 under the 
empirical reddening law described by Calzetti (2001). Figure 1 1 
shows the distribution of q values and the behavior of q with 
E(B — V). Interestingly, at large reddenings q never gets much 
above 1, suggesting that in these systems, the Lya emitting gas 
and the UV starlight are seeing the same amount of extinction. 
We do expect that galaxies with large Ly a optical depths will be 
censored out of our LAE sample. However, in dust-rich systems 
it appears that, if Lya escapes, it does so with very few resonant 
scatterings. This is consistent with models that involve strong 
winds, such as that proposed by Verhamme et al. (2008). On the 
other hand, at low reddenings, we see a large range of q values. 
Systems with q < 0 imply anisotropic emission, a top-heavy 
IMF, or a very young starburst, where the UV luminosity to SFR 
conversion breaks down. As expected, we find that the half-light 
radius and ^-factor are positively correlated, with a Spearman 
rank order coefficient of p = 0.35 (99.5% confidence). This 
correlation is shown in the right hand panel of Figure 12. A small 
size could lead to less homogeneity and thus more opportunities 
for Lya to undergo anisotropic radiative transfer. 

5. CONCLUSION 

Using broadband photometric data which extends from the 
rest-frame UV through to the near-IR, we have been able 
to measure the stellar masses, reddenings, and sizes for a 
sample of 63 luminous LAEs found in the HPS. Our fits 
demonstrate that, contrary to popular belief, Lya emitters are 
not exclusively low-mass objects. In fact, HPS-selected LAEs 
are quite heterogeneous, and are drawn from almost the entire 
stellar mass range of high-redshift galaxies. Moreover, there 
is a striking similarity between the mass function of LAEs 
and the mass function expected for the galactic star-forming 
population as a whole. This fact, and the lack of correlation 
between Lya luminosity and stellar mass, suggests that searches 
for Lya emission are excellent way of sampling a large fraction 
of the mass function of high-redshift star-forming galaxies. 

LAEs occupy a different part of stellar mass- SFR parameter 
space than that of galaxies found by other methods. Like the 
higher-mass submillimeter galaxies, LAEs fall above the main 
sequence of star-forming galaxies found by Daddi et al. (2007). 
This suggests that there is a different slope for the main sequence 
of star-bursting galaxies. Interestingly, LAEs do fall along the 


main sequence defined by Whitaker et al. (2012), though the 
~2 dex extrapolation required to reach their masses introduces 
significant uncertainty. Due to the various selection effects at 
work, the connection between the various classes of star-forming 
galaxies is murky at best. 

We also find that the range in observed ^-factors is depen- 
dent on the reddening, with the widest range of g- values oc- 
curring at low extinction. Interestingly, the observed values of 
q tend to unity as the reddening (or mass) increases, suggest- 
ing that in these objects, Lya photons are not undergoing a 
large number of scattering events. This strongly implies that 
winds are an important component in the making of high-mass 
LAEs. Furthermore, we find that the half-light radius and the 
^-factor are positively correlated, implying that Lya emission is 
enhanced in very small objects. 
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